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By means of ab initio calculations we study the effect of O-doping of Au chains containing a 
nanocontact represented by a Ni atom as a magnetic impurity. In contrast to pure Au chains, we 
find that with a minimun O-doping the 5 d xz , yz states of Au are pushed up, crossing the Fermi level. 

We also find that for certain O configurations, the Ni atom has two holes in the degenerate 3 d XZjyz 
orbitals, forming a spin S = 1 due to a large Hund interaction. The coupling between the 5 d xz , yz 
Au bands and the 3 d XZjyz of Ni states leads to a possible realization of a two-channel S = 1 Kondo 
effect. While this kind of Kondo effect is commonly found in bulk systems, it is rarely observed 
in low dimensions. The estimated Kondo scale of the system lies within the present achievable 
experimental resolution in transport measurements. Another possible scenario for certain atomic 
configurations is that one of the holes resides in a 3d z 2 orbital, leading to a two-stage Kondo effect, 
the second one with SU(4) symmetry. 
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I. INTRODUCTION 

In modern nanoscience, tailoring the electronic transport 
through atomic-size conductors has turned into a duty, 
as it is a powerful tool for detecting nanomagnetism. 
Atomic-size contacts can be experimentally obtained by 
different techniques, particularly in a mechanically con¬ 
trollable break junction (MCBJ) experiment, where the 
formation of one dimensional atomic chains of several 
kinds of elements is possible.® 5 Since the achievement 
of the first free-standing atomic chains of gold atoms in 
1998,®- the search for other elements that could also form 
atomic chains has been intense and active. For instance, 
in the last few years it has been possible to strengthen 
the bonds in a suspended chain and to achieve a higher 
probability of chain producibility by adding external ab- 
sorbates during the chain formation process. It is known 
that low-coordinated atoms are chemically more reactive 
than in bulk, 1 ^ thus, chains are expected to be even more 
reactive than nanoparticles, opening the possibility for 
molecular absorbates to dissociate, even at low tempera¬ 
tures. For instance, oxygen (O) atoms are expected to be 
incorporated in the chains, as predicted in several previ¬ 
ous works.®® 

Due to outstanding experimental achievements in the last 
decades, it is nowadays possible to design nanodevices as 
tools for detecting nanomagnetism by conductance mea¬ 
surements through atomic metal contacts. One can indi¬ 
rectly sense the presence of magnetism by detecting zero- 
bias anomalies, usually originated in the Kondo screening 
of the spin, if a magnetic impurity is bridging the contact. 
In previous works, Lucignano et al 13 and Miura et all 14 
studied the electronic structure and the Kondo conduc¬ 
tance through a Ni impurity embedded in a monoatomic 
Au wire. The Ni atom in the Au chain has two low-energy 


geometries: bridge (B) and substitutional (SUB). While 
in the B configuration, an usual one-channel Kondo 
(1CK) effect of spin S' = 1/2 was reported theoretically, 
no Kondo physics was predicted in the SUB configura¬ 
tion. This is due to the fact that in the SUB geometry 
(of higher total energy but probably accessible at large 
stress, as it is the case in a MCBJ experiment), the empty 
spin-down Ni state orbitals are a 3 d XZjyz degenerate pair 
with angular momentum projection \m\ = 1. In this ge¬ 
ometry the Ni impurity has S ~ 1 and the magnetic 
orbitals are ’spectators’ as they cannot be screened due 
to their orthogonality to the m = 0 Au conduction band. 
Indeed, transport experiments in Au chains indicate that 
the conduction channel is a single 6s band !® 16 

As it was previously mentioned, O impurities in the Au 
chain are expected during its real formation under an 
open atmosphere. A dramatic effect of the incorporated 
O impurities to the Au chains was found in our recent ab 
initio calculations. 17 They can modify the band struc¬ 
ture of the metallic host, pushing up the bands that 
are close to the Fermi level, thus establishing conduc¬ 
tion through the 5 d xz and 5 d yz electrons of Au. In these 
conditions, the spin of the Ni impurity, in the SUB ge¬ 
ometry, formed by the localized 3 d xz , yz electrons, can 
be screened by conduction electrons of the O-doped Au 
chain and Kondo physics is expected. A similar feature 
was obtained in O-doped Au chains containing a mag¬ 
netic S = 3/2 Co impurity!- 1 ^® Remarkably, the exam¬ 
ples of Ni-Au-0 and Co-Au-0 systems are expected to 
be realizations of Kondo models where an arbitrary spin 
S is screened by n channels, and were solved exactly in 
the 90’s by the Bethe ansatJ^Sl. For n = 2S the mod¬ 
els present Fermi liquids properties, while for n > 2S 
non-Fermi-liquid behavior is obtained and for n < 2S 
the systems are singular Fermi liquids. 21 The presence 
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of anisotropy render models not solvable by the B ethe 
ansatz and modify the low-temperature behavior . 18 19 22 
In this paper, we focus on a Ni magnetic impurity em¬ 
bedded, in a SUB configuration, within an O-doped Au 
chain. We study the effects of the presence of O impuri¬ 
ties on the symmetry of the Au conduction bands close 
to the Fermi level as well as on the spin state of the Ni 
impurity. This symmetry determines the nature of the 
screening of the impurity spin, giving rise to the possibil¬ 
ity of the system to exhibit some kind of Kondo effect. In 
particular, we found a stable configuration in which the 
model that describes the system corresponds to a spin 
5 = 1 screened by two conduction channels. Using pa¬ 
rameters coming from ab initio calculations, we build the 
model Hamiltonian which corresponds to a generalized 
Anderson impurity model. After a Schrieffer-Wolff trans¬ 
formation we prove that the effective low energy model 
is a S' = 1 , 2 -channel Kondo model. We analyze this 
model at zero temperature and give an estimation of the 
expected Kondo temperature (Tk) for this system. 

For other atomic configurations, corresponding to 
metastable states, the ab initio results suggest a six-fold 
degenerate ground state of the Ni atom, with one hole in 
a 3 d z 2 orbital and another one in either a 3 d xz or a 3 d yz 
orbital in an S = 1 state. This is analogous to the situ¬ 
ation of Fe iron(II) phtalocyanine molecules on Au(lll), 
where a two-stage Kondo effect takes place, screening 
first the spin 1/2 of the 3 d z 2 hole and then that of the 
orbitally degenerate remaining hole in an SU(4) Kondo 
effect 

The paper is organized as follows. In section [II] we pro¬ 
vide the details of our DFT first-principles calculations 
and show the band structure for the different hosts and 
the consequences of these structures on the systems with 
an embedded Ni impurity. In section IV we present the 
Kondo model, obtained by a Schrieffer-Wolff transfor¬ 
mation of the 5 = 1 two-channel Anderson model, and 
obtain the Kondo temperature. Finally, a summary and 
discussion are given in section |V| 


II. AB INITIO RESULTS: SYMMETRY OF THE 
Au CONDUCTION BANDS AND Ni HOLES 

We perform ab initio calculations based on density func¬ 
tional theory (DFT) using the full potential linearized 
augmented plane waves method, as implemented in the 
WIEN2K codeP^ The generalized gradient approxima¬ 
tion for the exchange and correlation potential in the 
parametrization of PBE 2 ^ and the augmented plane 
waves local orbital basis are used. The cutoff parameter 
which gives the number of plane waves in the interstitial 
region is taken as R mt * = 7, where K max is the 

value of the largest reciprocal lattice vector used in the 
plane waves expansion and R mt is the smallest muffin 
tin radius used. The number of k points in the Brillouin 
zone is enough, in each case, to obtain the desired energy 
and charge precisions, namely lO -4 Ry and 10 - 4 e, 



FIG. 1: (left) Electronic band structure and (right) or¬ 
bital decompositioned projected density of states, PDOS, of 
the Au monoatomic chain at the GGA equilibrium distance 
d & Au-Au ~ 4.9285 bohr. The Hubbard U parameter set in this 
calculations is U — 4eV. 


respectively. The muffin-tin radii were set to 2.23 bohr 
for Ni, 1.91 bohr for Au atoms and 1.69 bohr for the O 
impurities. In all the studied cases we consider a period¬ 
ically repeated hexagonal lattice with a = b = 15 bohr 
and the coordinate system is fixed in such a way that the 
chain axis is aligned with z. The a and b distances in the 
supercell were checked to be large enough to avoid artifi¬ 
cial interactions between the periodic replicas of the wire. 

As a first step, we analyze the case of a Au monoatomic 
chain and its band structure. Usually, DFT calculations 
yield a spurious magnetization in Au wires, due to 
self-interaction errors, which shift the 5 d XZjyz (\m\ = 1 ) 
bands of the Au wire up to the Fermi energy. One route 
to avoid this spurious result is to increase the Au-Au 
distance in the chains during the calculations . 13 14 We 
pick up another route to tackle this problem and, as 
suggested by Sclauzero and Dal Cor so ™ we include a 
Hubbard U correction of U = 4eV in the 5d-electron Au 
manifold. As it can be seen in Fig. [TJ the mentioned 
magnetization is completely suppressed in the chain at 
the equilibrium distance of d Au _ Au = 4.9285 bohr. 

In O-doped Au chains, we relax the Au-0 distance for 
the case of AuO diatomic chain (two-atom unit cell) 
and, afterwards, we take the same bond length, d Au _ 0 = 
3.625 bohr, for all the studied chains. We consider sev¬ 
eral O-dopings to find which the minimal amount of O 
is needed to push the projected 5 d XZjyz density of states 
of all Au atoms towards the Fermi level. In this way, 
the \m\ = 1 -symmetry conduction channel through all 
the Au atoms in the chain is opened. In all the studied 
cases - i.e.: 50% (AuO), 33.3% (A 112 O), 25% (AU 4 O) and 
14% (Au 6 0) - we include the same Hubbard U = 4eV 
in the 5d-electrons of Au and perform self-consistent cal- 
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FIG. 2: (Color online) (Top) 5 d xz , yz projected density of 
states of the Au atom located farthest away from the O atom. 
(Bottom) Schematic representation of the chains. The unit 
cell is shown for each case and the arrows indicate the Au 
atom whose projected density of states is shown in the top 
panel. 


dilations. We find that the 5 d XZjyz orbitals of Au, up 
to the third neighbors of the O impurity, cross the Fermi 
level, due to the large hybridization with the oxygen 2 p XrV 
states. This is shown in Fig.[2j where we plot the 5 d XZiyz 
density of states of the Au atom located farthest away 
from the O atoms, for some selected dopings. From this 
we can conclude that an O-doping of at least 14% is nec¬ 
essary to transform the Au bands into conducting chan¬ 
nels. We consider below a doping of « 19 % (3/16). 
Once the minimum amount of O present in the Au chains 
is set, we include the Ni impurity in our calculations. The 
question that arises then is, where to put the O atoms and 
the Ni impurity within the Au chain. We assume that 
the O atoms are not bonded to the Ni impurity and relax 
the Au-Ni bond length for a diatomic Ni-Au chain, giving 
dAu-m = 4.4805 bohr, in excellent agreement, within 
a 0.05%, with that obtained by Miura et al. £31. For a 
chain containig 12 Au atoms, 3 O atoms and 1 Ni impu¬ 
rity within the unit cell, we take the corresponding bond 
lengths as d Au _ Ni = 4.4805 bohr, d Au _ 0 = 3.625 bohr 
and d Au _ Au = 4.9285 bohr. In consequence, the Ni-Ni 
distance in the ^-direction is larger than 70 bohr and, 
therefore, we assume that the Ni-Ni magnetic interac¬ 
tion is negligible. 


To find out if the O’s position, with respect to the Ni 
atom, changes the spin-state of the Ni impurity or its 
hole’s symmetries, we explore several configurations, al¬ 
ways keeping the same amount of O in the Au chain. In 
tabled] we introduce the different cases studied and their 
corresponding Ni magnetic moments as well as the oc¬ 
cupation numbers (in electrons) obtained by integrating 
the minority-band separated in the different symmetries, 
within the Ni-mufhn-tin sphere. The last entry corre¬ 
sponds to a metastable solution, whose energy is 85 meV 
higher than that of the previous entry. Note that for the 
3d^-symmetry the maximal occupation could be as much 
as n(3d 2 z ) = 1, whereas n(3d xy , x 2_ y 2) = n(3d xz ^ yz ) < 2 as 
they both are 2-fold degenerate states. Based on the data 
shown in the table, we can distinguish between two ex¬ 
treme cases for the stable solutions: the one called sym 0- 
4 th -nn and the denominated sym 0-2 th -nn. In the former 
case, the Ni atom is placed in the center of the O-doped 
Au chain and it has two O atoms symmetrically located 
as fourth neighbors. The 3 d z 2 orbital is almost com¬ 
pletely filled and the two-fold degenerate orbitals 3 d XZjyz 
are almost empty. As the O atoms approach the Ni im¬ 
purity, the 3 d XZrVZ orbitals increase their filling while the 
3 d z 2 orbital moves towards the Fermi level decreasing its 
occupancy until it becomes almost half-filled in the case 
called sym 0-2 th -nn. In this last case, the Ni impurity 
has two O atoms at second neighbor positions. Note 
that for the metastable configuration, with one O atom 
as third and another as forth nearest-neighbors ( no-sym 
10-3 th -nn), the minority 3 d z 2 state is almost empty. The 
small energy difference (85 meV) between the metastable 
and stable configurations suggests a near degeneracy be¬ 
tween the triplet of the Ni 3d 8 configuration with two 
holes occupying the m = ±1 states (or 3 d xz , yz orbitals) 
and the two triplets with one hole in the m — 0 ( 3 d z 2 ) 
state and another one with m = ±1. In the calculations 
there is always a finite splitting, larger than 0.5 eV, pos¬ 
itive or negative of the peaks in the minority spectral 
density for m = 0 and m — ±1, but this is probably due 
to the exchange and correlation potential included in the 
DFT, since a partial occupation of one state pushes the 
other up in energy. 

In Fig. [3] we show the partial density of states projected 
onto the different symmetries of the Ni impurity, for two 
of the above mentioned cases. The Ni impurity develops 
a magnetic moment of (jlni = 1.37 fis in the sym 0~4 th ~ 
nn configuration and, as also found for the Ni impurity 
embedded in a Au monowirepl, the Au atoms at both 
sides of the Ni atom develop a small induced magnetic 
moment of around 0.1 8/j,b- The calculated Ni magnetic 
moment is fiNi = 1.33/i# in the sym 0-2 th -nn case, and 
the induced magnetic moment in the Au atoms located at 
both sides of the Ni impurity is around 0.31 /i#. The en¬ 
hancement of this last induced magnetic moment is due 
to the proximity of the O atoms (See sketch of Fig.[3jb)). 
In spite of the fact that the spin state of the Ni impurity 
in both configurations is almost the same, the hole’s sym¬ 
metries are different. It can be seen from Fig. |3ja) that 
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TABLE I: Symmetry-dependent d —band minority spin fillings of the Ni atoms (in electrons) for the selected studied cases. 
The color coding of the schematic representation of the chains is the one presented in Fig.[3|a). 


Case Unit cell’s schema n(3d z 2) n(3d xyx 2_ y 2) n(3d xz , yz ) 

symO-4 th -nn °G O 0*0 00*00 O c C C C 0.92 1.87 0.56 1.37 

sym 0-3 th -nn * 000©*0 0*0 0*0000 0.86 1.69 1.00 1.11 

sym 0-2 th -nn 0.63 1.60 1.21 1.33 

no-synn 10-2 th -nn *0060 0 * 0*0 00*000 0.67 1.58 1.18 1.31 

no-sym 10-3 th -nn * 0000*0 0*0 00*000 0.91 1.82 0.72 1.26 

idem above metastable ° 0000 * 00 * 000*000 0.27 1.53 1.67 1.26 


the empty spin-down Ni orbitals (holes) are the 3 d XZ:VZ , 
while the other orbitals remain occupied. The different 
band-fillings change, depending on the position of the 
O impurities. Upon approaching the O atoms towards 
the Ni impurity, the 3 d XZ:VZ orbitals begin to be filled, 
while the 3 d z 2 one begins to be unoccupied, as shown in 
Fig.§ b )- One also observes in this figure that the minor¬ 
ity 3 d z 2 state lies at the Fermi level and the density cor¬ 
responding to the 3 d XZr y Z states is split in two around the 
Fermi level, suggesting a near degeneracy of the m = 0 
and m = ±1 states as in the case when the O atoms are 
third and forth nearest-neighbors of Ni, already discussed 
above. In all the treated cases, the 3 d xy ^ x 2_ y 2 states are 
the most localized, as they lie perpendicular to the chain, 
thus having a small hybridization. Nevertheless, we ob¬ 
serve that the position of these filled degenerate orbitals 
might change depending on O proximity. 

In order to determine the anisotropy constant D , taken 
into account in the term Hp = DMl of the Hamilto¬ 
nian defined in Eq. [l] of Section |IV[ we introduce the 
spin orbit interaction in the Ni impurity and perform 
self consistent calculations within the fully relativistic 
approximation. For the sym 0~4 th ~nn case we obtain 
a magnetocrystalline anisotropy energy, MAE, equal to 
E\\ — Ersj 10 meV, where E\\ (E_ l) is the total energy 
corresponding to the total magnetization pointing along 
(perpendicular to) the chain’s axis. Therefore, we ob¬ 
tain that the easy magnetization axis is perpendicular 
to the chain. To calculate the value of D we proceed 
as in Ref. EH The expectation value of the anisotropy 
in the state 111) of maximum projection in the 2 direc¬ 
tion is clearly (11\Hd\11) = D. By rotating this state, 
one obtains that the state of maximum projection in 
the x direction is \llx) = |10)/\/2 + (111) + |1 — l))/2, 
and (11x\Hd\11x) = D/2. From the energy difference 
D ~ 20 meV. 


III. EFFECTS OF CORRELATIONS INSIDE 
THE Ni 

In this section we investigate the effects of correlations 
in an isolated 3d 8 configuration on the relative stability 
of the two possible ground states discussed above: spin 
triplet with total angular momentum projection L z = 0 
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(a)Case sym 0~4 th -nn : Two O atoms are 
symmetrically positioned as Ni-fourth neighbors. 
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(b)Case sym 10-2 th -nn : Two O atoms as 
Ni-second neighbors. 

FIG. 3: Ni impurity in a 19% O-doped Au chain. In each 
case, the unit cell and the density of states projected onto 
different symmetries are shown. 


with one hole with m = 1 and the other with m = — 1, 
or two triplets with L z = ±1 (one hole with m = 0 and 
the other with m = ±1). We also provide an alterna¬ 
tive estimation of the anisotropy parameter D. Since the 
GGA underestimates correlations that affect the orbital 
polarization of the d statesthis calculation is an 
important complement to the ab initio results presented 
above. 
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We have solved exactly the 45 x 45 matrix of the Hamil¬ 
tonian corresponding to the 3d 8 configuration, contain¬ 
ing crystal-field splitting and all correlations inside the d 
shell, as described for example in Refs. 1311321 In a later 
step, to calculate D , we include the spin-orbit coupling 
Hsoc = A JAl* • Si. We have considered the Coulomb 
integrals F 2 = 0.16 eV, F 4 = 0.011 eV, and A = 0.08 
eV from a fit of the low energy spectra of late transition 
metal atoms. 

For the crystal-field splitting, we took two different sets of 
parameters. From the ab initio results presented in Fig. 
3 (a) for the minority states, we estimate -1.85, -1 and - 
0.3 eV for the on-site energy of the states with m = ±2, 0 
and ±1 respectively. For the latter, which is unoccupied 
we have subtracted a correlation term included in the 
GGA, which we estimate to be 0.5 eV. In the second 
set we assume that the m = 0 and m = ±1 states are 
degenerate, with the same average distance to the m = 
±2 states as before. This is equivalent to take -1.2 eV for 
the energy of the electrons with x 2 — y 2 and xy symmetry 
and 0 for the remaining three. Note that in the last case, 
in the absence of interactions, the above mentioned two- 
hole possible ground states with L z = 0 and L z = ± 1 are 
degenerate. We find that interactions stabilize the triplet 
with L z — 0 by 0.30 eV. For the first set of parameters, 
the stabilization energy for this triplet, which is 0.7 eV 
in the absence of interactions, is increased to 0.88 eV in 
the present calculation. The difference 0.18 eV can be 
ascribed to the effect of correlations. 

When spin-orbit is included we obtain D = 8.5 (16.7) 
meV for the first (second) set of parameters. The sign 
and order of magnitude agree with the ab initio results. 


IV. Ni IMPURITY WITHIN AN O-DOPED Au 
CHAIN: KONDO EFFECT 

According to the results of the previous section, different 
model Hamiltonians can be formulated depending on the 
O position relative to the Ni impurity. In general, for 
the Ni 3d 8 configuration, the triplet ground state with 
L z = 0 is the most stable in the ab initio calculations 
and is further stabilized by correlations. Therefore in 
this section, we analyze this case in more detail. The 
situation with two degenerate triplets with L z = ±1 is 
very similar to the physics of an organic molecule in a 
metallic substrate and will be addressed to in Section El 
To be specific, the system we are interested in, repre¬ 
sented by the first row of table [l| namely sym 0-J^ h -nn^ 
is the one in which the symmetries of the Ni holes are 
well defined in the orbitals a = yz,xz. Moreover, this 
case presents a larger hybridization with the host as it 
can be seen in the broader density of states of these or¬ 
bitals, see Fig. [3j leading to a larger Kondo scale. Tak¬ 
ing into account the ab initio results presented in the 
previous section, the charge fluctuations in the Ni impu¬ 
rity are mainly due to the interchange between the state 
with two holes and virtual excitations to states with only 


one hole in either the 3 d yz or 3 d xz orbital (See Fig. [4|. 
Furthermore, both configurations are magnetic with spin 
5 = 1 and S = 1/2 for two and one holes respectively. 
The spin-orbit coupling (SOC) in the Ni atom induces a 
splitting D between the projections 5^ = 1 and 5^=0 
of the triplet that belongs to the total spin 5 = 1. 



FIG. 4: Schematic representation of Ni-holes and their tun¬ 
neling to the Au-leads for the case sym 0-Jf h -nn. 

With these information, we propose a Hamiltonian that 
describes the system, and that is given bj j33] 


H = ^(E 2 + DM%)\M 2 ){M 2 \ + E a \aM 1 )(aM 1 \ + 

M 2 olMx 

y! e vk c ];koi(j C vkacr + (1) 

vkOLCJ 

E E V I/ (lM 2 \^t Ml a)(\M 2 )(aM 1 \c, kai 7 +n.c.), 

Mi M 2 auka 

where Ei and Mi represent the energies and the spin 
projections along the chain, chosen as the quantization 
axis, of states with i — 1,2 holes in the 3d shell of the 
Ni impurity. The state with two holes and maximum 
spin projection is denoted by |1) = d} xz ^d) yz ^ |0), where 
|0) represents the 3d 10 configuration and the operator 
F ao . creates a hole with symmetry a = xz,yz and spin 
projection a. The states with one hole in the Ni atom 
can be constructed by removing an a hole. The other 
relevant states with two and one holes can be obtained 
by using the spin lowering operators. 

The operator c^ kaa creates a hole in the 5 d shell of the 
Au atom with symmetry cq where v = L,R denotes the 
left or the right side of the Ni atom, respectively. The 
hopping V v characterizes the tunneling between the Ni 
and Au states. 


Tine Kondo regime. 

To describe the spin fluctuations, ie, the Kondo regime 
in which the charge fluctuations are frozen, we perform 
a Schrieffer-Wolff transformation. The effective low en¬ 
ergies’ Hamiltonian is then 
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Heff — 6 ukC tkaa CukacT + 

uau 

^ |U| 2 f |l)(l| + |-l)(-l| |o)(o| ^ , 

^ 2 1 El-D-E 2 + El - E 2 ' cm/ + 

Otis 

J w E + S v s i ) (2) 

a a 

where h cav = Y, ka c[ ka<J c vkoia and Jy = El ^_ D , J± = 

!/2 (eM^d + AJ and l^l 2 = {\Vl\*+\Vr\ 2 )/2. The 
middle term represents a renormalization of the local en¬ 
ergies and the potential scattering and are unimportant 
for understanding the Kondo effect. The final term rep¬ 
resents the exchange interaction between the impurity 
spin 5 and the corresponding ones of the O-doped Au 
bands s a , being Jy and J± the longitudinal and perpen¬ 
dicular couplings. This exchange term corresponds to the 
anisotropic Kondo Hamiltonian, Hk 

Hk = J„E^4 + ^E(« + «) (3) 

ot a 

Notice that in the absence of any anisotropy interaction, 

D = 0, the model reduces to J 5 • ( s yz + s xz ), which 
is a spin-1 two-channel, n — 2, Kondo model with J = 
e^-e 2 > * s a Fermi liquid pSEll 

The Kondo temperature for our model can be estimated 
from the usual expression Tk ~ e i/2jp c (E F ) - m w pi c ] 1 
p c (Ep ) corresponds to the value of the density of states 
of the conduction band at the Fermi level, Ep. The 
only change with those associated to the spin-1/2 one- 
channel model enters in the change of the exchange cou¬ 
pling J, a half in our case as compared to the spin-1/2 
one-channel.^ 

For a numerical evaluation in terms of the ab initio pa¬ 
rameters we use the fact that Jp c (Ep) = being 

A = ttV 2 p c (Ep) and = E\ — E^ ~ —E xz . We es¬ 
timate A from the half maximum width of the 3d a den¬ 
sity of states, which corresponds to 2A ^ 0.23 eU, and 
\E XZ \ ~ 0.23 eV from its energy position. Thus, we es¬ 
timate a Kondo temperature close to T K ~ 400A, and 
therefore the zero-bias anomaly should be observed in 
transport measurements. 

As stated in Section [ilj we find that the value of the 
anisotropy constant D is of the order of D ~ 15 meU ^ 
180 K. 

Studies of similar models by the numerical renormaliza¬ 
tion group suggest that for D < Tk (as in this case) an 
energy scale D* much smaller than D is dynamically gen¬ 
erated. For example, for the underscreened 5 = 1 Kondo 
model with positive D, D * « exp[—cfTV/Z)) 1 / 2 ]. While 
for D = 0 the conductance reaches its maximum value 
at T = 0, when D ^ 0 and for T < D* the conductance 
drops abruptly? 22 A similar D* was found in an under¬ 
screened 5 = 3/2 model which displays non-Fermi liquid 
behavior for T < D* . 18 


Thus, while for temperatures above D*, the model ex¬ 
hibits the features of the isotropic 5 = 1 two-channel 
model, and the splitting of the triplet states can be in 
practice neglected for T D *, at temperatures be¬ 
low D* one expects a change of regime and a drop in 
the conductance. This is a general feature of the mod¬ 
els for which the spin degeneracy of the impurity is re¬ 
moved,^ as in the underscreened 5 = 1 Kondo model 
with anisotropy. 22 


V. SUMMARY AND DISCUSSION 

By means of ab initio calculations we determined that 
the minimal amount of O-doping needed in a Au chain 
to push the 5d-bands towards the Fermi level is 14%. 
With this amount of O atoms, the 5 d xz ^ yz orbitals of all 
the Au atoms present in the chain are conduction bands, 
normally below the Fermi level when not doped. Within 
this minimal concentration of dopants, we studied the 
effect of the O atoms on the spin state of a Ni mag¬ 
netic impurity embedded in the O-doped Au chain, as 
well as the symmetry properties of the unoccupied bands 
(holes). We showed that when the O-atoms are fourth 
neighbors of the Ni impurity, the 5 = 1 spin state of the 
Ni atom comes from its 3 d XZjyz empty states. This is a 
triplet with total angular momentum projection L z = 0. 
Upon approaching the O atoms towards the Ni impu¬ 
rity, the 3 d xz , yz begin to be filled, while the 3d 2 starts 
to be unoccupied, preserving the 5 = 1 spin. Therefore 
a transition to a Ni 3d 8 configuration with two triplets 
with L z = ± 1 becomes possible, although correlations 
stabilize the L z = 0 triplet. 

After the description of the system and of the first prin¬ 
ciples calculations, we introduced a mixed valence model 
for the Ni-O-Au chain in the substitutional configura¬ 
tion for Ni, assuming that the Ni impurity for two holes 
is in the L z = 0 state. This can be viewed as a gen¬ 
eralized Anderson impurity model, which includes not 
only the charge fluctuations (between 3d 8 and 3d 9 Ni 
configurations) but also the spin ones. By means of a 
Schrieffer-Wolff transformation we prove that the model 
can be mapped onto a two-channel spin-1 Kondo model. 
The SU(2) symmetry of the conduction channels of the 
O-doped Au chain allows for a full screening of the spin 
of the Ni impurity. We found that the associated Kondo 
temperature is experimentally accessible and that it is 
of the order of Tk ~ 400 K. Therefore, at low enough 
temperatures as compared with the Kondo one, the sys¬ 
tem behaves as a Fermi liquid and a zero-bias anomaly 
should appear in transport measurements. We also found 
that even if the low energy model corresponds to an 
anisotropic one, the expected behavior should be similar 
to the isotropic case due to the fact that the anisotropy 
constant is found to be smaller than the isotropic Kondo 
temperature. A small energy scale D* is dynamically 
generated, below which the conductivity should drop. 
We hope that this work stimulates an accurate many- 
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body calculation of D* which might be compared with 
the D* & exp [—c(Tjf /D) 1 11 / 2 ] behavior obtained in similar 
modelsP^ST 

As pointed out in the seminal work by Nozieres and 
Blandinp^ if the relation n = 2S is satisfied (as it is 
actually the case), then, at temperatures below the char¬ 
acteristic one given by the Kondo temperature, T <C TV, 
the physics of the system is expected to be the same that 
for a Fermi liquid. This is an example of a fully screened 
higher spin Kondo effect. Most of the transition metals 
in bulk conducting materials are described by this kind 
of totally screened Kondo phenomena.^ 5 However, this is 
not the usual behavior in low dimension and therefore 
our model represents a possible experimental realization 
of this kind of physics. Other low dimensional exam¬ 
ples of a two-channel spin-1 Kondo effect can be found 
in quantum dots (QD). Specifically, Pustilnik and Glaz- 
man 37 !! s t u qi e d the possibility of having a spin-1 QD 
coupled to two screening channels which fully screen the 
QD spin. However, they found that even with a small dif¬ 
ference between the two antiferromagnetic coupling con¬ 
stants (which corresponds to a real situation), the physics 
of the QD is expected to be dominated by the under¬ 
screened one-channel spin-1 Kondo model, leading to a 
two-stage Kondo effect Fortunately, our model based 
on the O-doped Au chain has the advantage of having a 
protected SU(2) symmetry between the two conducting 
channels due to the cylindric symmetry. 

The case with two degenerate triplets in the 3d 8 con¬ 
figuration, in which one hole occupies the 3d orbital of 


z 2 symmetry and the other one can be occupied by an¬ 
other hole with xz or yz symmetry is completely analo¬ 
gous to the case of iron(II) phtalocyanine molecules on 
Au(lll). 23 In this case, due to the different hybridization 
of the Fe orbitals with the substrate (as in our case, the 
hopping is larger for the z 2 orbital), a two-stage Kondo 
effect takes place. First, the spin of the z 2 orbital is 
screened. It is difficult for us to estimate the Kondo scale 
for this screening due to the uncertainty in the position 
of the 3 d z 2 state. The remaining degrees of freedom, 
the orbital (xz or yz) degeneracy and spin 1/2 of the re¬ 
maining hole are screened in a rather exotic SU(4) Kondo 
effect. The effects of the orbital degeneracy become evi¬ 
dent when the molecules are assembled in a lattice SSEU 
It would be certainly interesting to find an SU(4) Kondo 
effect in nanoscopic chains. 
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